Long-range entanglement in the XXZ Heisenberg spin chain after 

a local quench 

Jie Ren^* and Shiqun Zhu^^ 
^Department of Physics and Jiangsu Laboratory of Advanced Functional Materials, 
Changshu Institute of Technology, Changshu, Jiangsu 215500, China 
"^School of Physical Science and Technology, 
Suzhou University, Suzhou, Jiangsu 215006, China 
(Dated: June 2, 2010) 

Abstract 

The long-range entanglement dynamics of an one-dimensional spin-1/2 anisotropic XXZ model 
are studied using the method of the adaptive time-dependent density-matrix renormalization-group. 
The long-range entanglement can be generated when a local quench on one of boundary bonds is 
performed in the system. The anisotropic interaction has a strong influence both on the maximal 
value of long-range entanglement and the time of reaching the maximum long-range entanglement. 
The local coupling has a notable impact on the long-range entanglement, but it can be neglected 
in the time of reaching the maximal long-range entanglement. 
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I. INTRODUCTION 



Entanglement generation and distribution is one of the important problems in performing 
mantum-information tasks, such as quantum computation and quantum teleportation [1- 
|. Many results showed that entanglement existed naturally in the spin chain when the 



temperature is at zero 



5[. It is discouraging that the entanglement in many spin systems 
is typically very short ranged. It exists only in the nearest neighbors and the next nearest 
neighbors js]. It is interesting that some schemes for long-range entanglement, such as 
exploiting weak couplings of two distant spins to a spin chain were proposed UB- These 
methods have limited thermal stability or very long time scale of entanglement generation. 
In recent years, many researches has shifted the focus of the study in the dynamics of 



entanglement 



10Nl4|. The dynamics of long range entanglement are obtained by a time 



quench of magnetic field 15|, |l6| . The end-to-end entanglement can be shared across a chain 



of arbitrary range 



lOj. It is excited that the most common situation studies so far concerns 



17 



19|. In Ref. (l9|, it showed 



a sudden quench of some coupling of the model Hamiltonian 
that the long-range entanglement can be engineered by a non-perturbative quenching of a 
single bond in a Kondo spin chain with impurity. This is the first example that a minimal 
local action on a spin chain can generate long-range entanglement dynamically. It would be 
interesting to investigate the possibility of producing long-range entanglement in the XXZ 
Heisenberg spin chain without impurity. 

It is well-known that the Hamiltonian of an opened chain of N spin-1/2 system with 
nearest-neighbor XXZ interaction is given by 



N-l 



(1) 



i=l 



where S°'{a = x, y, z) are spin operators on the z-th site, is the length of the spin chain. 
The parameter A denotes the couplings in the z-axis.This model can be realized in the 



Josephson-j unction 



20] and optical lattices 2l|, |22| 



In the paper, the long-range entanglement in the XXZ Heisenberg spin chain after a 
local quench is investigated. An opened boundary condition (OBC) is assumed because 
the antiferromagnetic Heisenberg spin chain with OBC can be achieved artificially in the 
experiment 23|, and the coupling J = 1 is considered for simplicity. In section II, the local 



quench is presented. The concurrence is used as a measurement of the entanglement. In 



section III, a more general situation and a single quench and its effect on the end-to-end 
qubits entanglement are analyzed. The effects of the anisotropic interaction and system size 
on the end-to-end qubits entanglement is also studied. The robustness of the entanglement 
against an increase in temperature is investigated. A discussion concludes the paper. 

II. LOCAL QUENCH AND ENTANGLEMENT MEASURE 

In the paper, the system is assumed to be in the ground state Gsq of Hq initially. A local 
quench change of the coupling between the first qubit and the second qubit with the same 
anisotropy interaction A in Eq. (1) is performed. The Hamiltonian of the system modifies 
to 

N-l 

H, = MS^.SI + SfS', + ASfSI] + J[S:s:^, + SfSl, + ASISI^,]. (2) 

i=2 

Since [Hq, H]\ ^ 0, the the ground state Gsq of Hq is not one of the eigenstates of Hi. The 
state of the system will evolve as 

= exp-^-^^* Gsq. (3) 

In the paper, the concurrence is chosen as a measurement of the pairwise entanglement 
[3I. The concurrence C is defined as 

Ci^N = max{Ai - A2 - A3 - A4, 0}, (4) 

where the quantities Aj(i = 1, 2, 3, 4) are the square roots of the eigenvalues of the operator 
Q = Puio-f (g) o-^)Pi N^'^i ® ^n)- They are in descending order. The case of Ci^n = 1 
corresponds to the maximum entanglement between the two qubits, while Ci^n = means 
that there is no entanglement between the two qubits. 

III. LONG-RANGE ENTANGLEMENT DYNAMICS 

It is known that it is hard to calculate the dynamics of entanglement because of the lack 
of knowledge of eigenvalues and eigenvectors of the Hamiltonian. For models that are not 
exactly solvable, most of researchers resort to exact diagonalization to obtain the ground 



3 



state for small system size. The method is difficult to be applied for large system size 
> 20. For a large system, the adaptive time-dependent density- matrix renormalization- 
group can be applied with a second order Trotter expansion of the Hamiltonian as described 



m 



24 



25| . In order to check the accuracy of the results of the adaptive time-dependent 
density-matrix renormalization-group, the results of exact diagonalization can be considered 
benchmark for a small size system. 
The numerical error of the adaptive time- dependent density-matrix renormalization- 
group comes from the discarded weight and the Trotter decomposition. The error of dis- 
carded weight is dependent on the data precision and truncated Hilbert space. The error 
of the Trotter decomposition is relied on Trotter slicing. In our numerical simulations a 
Trotter slicing 5t = 0.05 and Matlab codes of the adaptive time- dependent density-matrix 
renormalization-group with double precision are performed with a truncated Hilbert space 
of m = 100. In turns out that a typically discarded weight of 6p < 10^^. The error of 
Trotter decomposition 6 oc These can keep the relative error 6C in C below 10^'^ for 

a chain of = 60 sites with time t < 50/ J. 

In the adaptive time-dependent density-matrix renormalization-group, it is hard to calcu- 
late the reduced density matrix pi^N- By making use of the relation between the correlation 
and the reduced density matrix, the reduced density matrix can be expressed as 

Pl,N = ^[h,N+ Yl (5) 

where cr^(a = x,y,z) are Pauli operators and Ii^n is identity matrix. 

We obtain the reduced density matrix pi jv by Eq. (5), then calculate the entanglement 
between two-ends qubits. The end-to-end entanglement Ci^^ is plotted as a function of 
the anisotropic interaction A and time t with Ji = —0.01 in Fig. 1(a). It is seen that 
the long-range entanglement Ci^n can be generated after a short time. When the time t 
increases, the entanglement reaches the maximal value and then disappears quickly. That 
is, there is a peak in the entanglement. The peak of the end-to-end entanglement is plotted 
as a function of the anisotropic interaction A in Fig. 1(b). It is seen that the peak in 
long-range end-to-end entanglement increases with the anisotropic interaction A. It reaches 



the maximal value when the anisotropic interaction A = 1 j26|, 1271]. With the anisotropic 



interaction A increases further, the height of the peak decreases. 
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The relation between the time of reaching the maximal end-to-end entanglement labeled 
by Tjnax and the anisotropic interaction A can be seen in the inset of Fig. 1(b). The time of 
reaching the maximal long-range entanglement decreases when the anisotropic interaction 



A increases, while they are also not a linear relationship. Similar to Refs. [10|, Il9|], the 
system generates end-to-end entanglement periodically. In our simulations, the error of the 
simulations increases when the time increases. Since the coherent time of the system is not 
very long, the long-distance entanglement is not shown when the time t>50/J. It is shown 
that a sharp drop of long-distance entanglement occurs at A = —0.5. Furthermore, the 
time of the long-range entanglement reaching the maximal value goes up drastically. This 
is similar to the anomalous behavior appeared in Ref. 27 1. 

The end-to-end entanglement Ci^n is plotted as a function of the interaction Ji and the 
time t with A = 1 in Fig. (2). The size of the system = 20 is chosen. It does not 
include the case of Ji = 0. The long-range entanglement Ci a? can be generated after a 
short time. There is a peak in Ci^n when Ci^n is plotted as a function of time t. It is seen 
that the influence of the coupling Ji on the maximal long-range entanglement is relatively 
large. The peak of the end-to-end entanglement increases when the coupling interaction Ji 
increases, and reaches the maximal value when Ji = —0.1. It seems that the changing of local 
interaction from antiferromagnetic to ferromagnetic may enhance the entanglement creation 
When Ji increases further to Ji > 0.32, the long-range entanglement disappears. It is 
interesting that the coupling interaction Ji has a relatively small impact on the time when 
the end-to-end entanglement generates and reaches the maximal value. 

It is easy to obtain that [HQ,J2iLi ^i] = [Hi,J2iLi Si] = 0. It means that the ground 
states of Hq is a total singlet Stot = 0. The Hamiltonian Hi has invariant on every excitation 
subspaces. During the evolution, the boundary spin at z = 2 will have a strong tendency 
to form a singlet pair with its only nearest neighbor i = 3 on the right-hand side. This is 

similar for spin pairs (4, 5) and (6, 7), , etc. Then the two-end spin qubits form a singlet. 

The local interaction Ji decides the ability of forming singlet for even bands. Thus, the 
long-distance entanglement creates. 

The thermalization and relaxation during the period of generating long-distance entan- 
glement are neglected because the dynamical time scale is quite short. When initial state is 
taken to be the relevant thermal state, the end-to-end entanglement is plotted as a function 
of temperature in Fig. 3 after bond quenching for different system sizes. It is shown that 
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the entanglement vanishes when kT > 1.16 with the size of = 8, and kT > 1.06 with the 
size of = 10. It seems that our scheme is quite robust against temperature and is similar 
to that in Ref. jl9|. The initial thermal state does not have long-distance entanglement in 
Eq. (1), but the energy gap between the ground state and the first and other excited states 
are larger than the system with impurity. This leads to our scheme is quite robust against 
temperature. 

It is interesting to investigate the long-range entanglement creation even for very long 
chain of large size. The maximal value of the end-to-end entanglement Ci^n is plotted 
as a function of the size of the system for different anisotropic interaction A and different 
interaction Ji in Fig. 4(a). It is found that the maximal value of the end-to-end entanglement 
Ci^jv decreases when the size of system increases. When the anisotropic interaction A = 1 
and the interaction Ji = —0.1, the maximal value is 0.5481 for = 40 and 0.4547 for 

= 60. It is shown that the maximal value of the end-to-end entanglement label by ^(A^) 
is given by 10 1 



^(A^) ~ 1.35N-^/\ 



(6) 



In Ref. [l0|, ^(40) ~ 0.3947 for A^ = 40, and ^(40) ~ 0.3448 for A^ = 60. While in the 
Kondo regime, ^(40) can be as large as 0.7 for A^ = 40. It seems that our results are 
also quite good. The long distance entanglement can be obtained by performing a local 
quench in either one of two boundary bands. The time of reaching the maximal end-to- 
end entanglement labeled by Tmax is plotted as a function of the size A^ of the system for 
different anisotropic interaction A and different interaction Ji in Fig. 4(b). The time of 
the end-to-end entanglement reaching the maximal value Tmax is linear increase when the 



size A^ increases 
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27] . It is seen that the slope of the line is almost dependent on the 



anisotropic interaction A. 

In the paper, adiabatic quenches are studied and the decoherence effects in the system 
are ignored. It is noted that the XXZ spin chains can be realized in the Josephson-junction 
array 20|and optical lattices 2l|, |22|. In Josephson-junction array, the interactions and 



the anisotropy can be modulated by varying voltages [28|]. The Josephson-junction array 
system does not have significant decoherence over our time-scales {Tmax — N/J) [29|. The 



anisotropic interaction quenches can be achieved successfully with long decoherence time in 



the optical lattices experiment 



22| 



IV. DISCUSSION 



By using the method of the adaptive time-dependent density-matrix renormahzation- 
group, the time evolution of the entanglement in a one-dimensional spin-1/2 anisotropic 
XXZ model is investigated when a local quench is performed in the system. The local 
quench is a abrupt change of interaction between the first qubit and the second qubit. 
The dynamics of pairwisc entanglement between the two ends qubits in the spin chain is 
studied. The entanglement of the two-ends spin qubits can be created after the local quench 
is performed.The time when the long-range entanglement generates decreases in a nonlinear 
manner with the anisotropic interaction increases. It reaches the maximal value when the 
anisotropic interaction A = 1. The maximal value of the end-to-end entanglement increases 
with the coupling interaction Ji increases, and reaches the maximal when Ji = —0.1. It 
is interesting that the coupling interaction has a relatively small impact on the time when 
the long-range entanglement reaches the maximal value. This phenomenon may be used 
to control the dynamics of the entanglement by varying the anisotropic interaction and the 
local quench interaction of the Heisenberg spin chain. 
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FIG. 1: (a). The long-range entanglement Ci^n is plotted as a function of anisotropic interaction 
A and time t with Ji = —0.1. The size of the system is A'^ = 20. (b). The maximal value of 

the cnd-to-cnd entanglement is plotted as a function of the anisotropic interaction A. The inset 
shows the time of reaching the maximal end-to-end entanglement labeled by Tjnax as a function of 
anisotropic interaction A. 
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FIG. 2: The long-range entanglement Ci^n is plotted as a function of interaction Ji and time 
t when there is local quench of the interaction Ji. The size of the system is = 20 and the 
anisotropic interaction A = 1. It does not include the case of Ji = 0. 



10 



1.0- 




0.0 0.3 0.6 0.9 1.2 1.5 

kT 



FIG. 3: The maximal value of end-to-end entanglement is plotted as a function of temperature 
after band quenching for different sizes. 
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FIG. 4: (a). The maximal value of the long-range entanglement Ci^n is plotted as a function of 
the size of the system for different anisotropic interaction A and different interaction Ji . (b) . The 
time of reaching the maximal Ion g-range cntcinglGrnent labGlcd. by T^nax 

is plotted function of 

the size of the system for different anisotropic interaction A and different interaction Ji. The red 
lines are fixed lines. 
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